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Abstract. - The time-dependent scaling of the two-time autocorrelation function of spin 
systems without disorder undergoing phase-ordering kinetics is considered. Its form is shown 
to be determined by an extension of dynamical scaling to a local scale-invariance which turns 
out to be a new version of conformal invariance. The predicted autocorrelator is in agreement 
with Monte-Carlo data on the autocorrelation function of the 2D kinetic Ising model with 
Glauber dynamics quenched to a temperature below criticality. 



Understanding the kinetics of phase-ordering after a rapid quench from an initial disordered 
state into the ordered phase has since a long time posed a continuing challenge (see [1-4] for 
reviews). A key insight has been the observation that many of the apparently erratic and 
history-dependent properties of such systems can be organized in terms of a simple scaling 
picture [5]. This means that there is a single time-dependent length-scale L(t) which is 
identified with the typical linear size of ordered clusters. It turns out that the ageing behaviour 
is more fully revealed in observables such as the two-time autocorrelation function C(t,s) or 
the two-time linear autoresponse function R(t, s) defined as 



C(t,s):=(<j>(t)<f>(s)) , R(t,s)- ! 



5h(s) 



(1) 

h=0 



where <p(t) denotes the time-dependent order-parameter, h(s) is the time-dependent conjugate 
magnetic field, t is referred to as observation time and s as waiting time. One says that the 
system undergoes ageing if C or R depend on both t and s and not merely on the difference 
r = t — s. These two-time functions are expected to show dynamical scaling in the ageing 
regime t,s ^> t m icro and t — s ^> £ m i cro , where t mlcIO is some microscopic time scale. Then 

C(t, s) = M^f c (t/s) , R(t, s) = s- l - a f R {t/s) (2) 

such that the scaling functions fc,n{y) satisfy the following asymptotic behaviour 

fc(y) - y- Xc/z , f R (y) - 2T Wz (3) 
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as y — y oo and where Ac and A/j, respectively, are known as the autocorrelation [6,7] and 
autoresponse exponents [8] and z is the dynamical exponent, defined through L(t) ~ i 1 / 2 . 
Throughout, we consider simple ferromagnets without disorder and with a non-conserved 
order-parameter. Then z = 2 is known [9]. For spin systems with short-ranged equilibrium 
correlators (e.g. the d > 1 Glauber-Ising model) it has been checked in detail that a — 1/z = 
1/2 [10,11], but the closed-form and ad hoc OJK approximation gives a — (d— 1) /2 [12,13]. The 
exponents Xc.r are independent of the equilibrium exponents and of z [2,3, 14]. Although 
the equality Ac = A^ had been taken for granted (reconfirmed in a recent second-order 
perturbative analysis of the time-dependent Ginzburg-Landau equation [15]), counterexamples 
exist for long-ranged initial correlations in ageing ferromagnets [8] and in the random-phase 
sine-Gordon model [16]. For short- ranged initial correlations, dynamical scaling together with 
Galilei-invariance at temperature T = are sufficient for Ac = A# [17]. 

We are interested in the form of the scaling functions fc,n(y)- Indeed, it is known that 
for any given value of z there exist infinitesimal local scale-transformations t i— ► (1 + e) z t, 
rH (1| e)r with an infinitesimal e = e(t, r) which may depend on both time and space [18]. 
Furthermore, the local scale-transformations so constructed act as dynamical symmetries of 
certain linear field equations which might be viewed as some effective renormalized equation of 
motion. From the assumption that the response functions of the theory transform covariantly 
under local scale-transformations, the exact form of the scaling function fn(y) is found [17-19] 

f R {y) = r o2 / 1+a '- A -/ z (y ~ l)- 1 -*' (4) 

where a' is a new exponent [17] and ro is a normalization constant. Eq. J3J with a = a' is 
recovered in many spin systems quenched to a temperature T < T c and whose dynamics is 
described by a master equation [4,17-21], but in models such as the II? Glauber-Ising model 
at T = [17] or the OJK approximation of phase-ordering [12, 13, 15] eq. Q holds but with 
a ^ a! . If a phase-ordering system is also Galilei- invariant at T — 0, then fn(y) is independent 
of both the thermal and the initial noises [17]. 

While the scaling form of the autoresponse function thus seems to be understood, the 
problem of finding the scaling function fc(y) of the autocorrelation function appears to be 
considerably more difficult. A by now classical attempt recognizes that for T < T c , temper- 
ature should be irrelevant [1] and hence sets T = 0. Building on the Ohta-Jasnow-Kawasaki 
approximation (see [1]) in the kinetic 0(n)-model one introduces an auxiliary field for which 
a gaussian closure procedure is assumed. This leads to [22-24] 

, . n r (1 n + l\l 2 / 4y \ d/i (l 1 n + 2 / Ay \ d/2 \ . . 

fc,BMy) = ^[B^—)\ [ty-Tir) ^[rr—'lty-TW) ) (5) 

where B is Euler's beta function and %F i a hypergeometric function. However, this closed form 
implies Ac = d/2 which only holds in certain limiting cases (for example, Ac = d/2 + an -1 , 
to leading order in n, in the 0(n)-model and with a known value of a > [25]. See also [26]). 

Here we investigate to what extent fc{y) m ay be determined from a local scale-invariance 
(LSI). We concentrate on phase-ordering where T < T c and thus z = 2 [9]. The group of 
local scale-transformations is then the Schrodinger group [27, 28] which for example arises as 
the maximal kinematic group of the free Schrodinger (or diffusion) equation. In particular, 
the Schrodinger group contains dilatations with z — 2 and Galilei-transformations. For local 
theories, there is a Ward identity such that these two symmetries imply full Schrodinger- 
invariance [29]. However, Galilei-invariance is incompatible with thermal or initial noises. 
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We consider a coarse-grained order-parameter cj>(t, r) satisfying a Langevin equation [1,30] 

= - D 5 2L Dv{tm r) + ??(t , r) (6 ) 

where TL is the classical Hamiltonian, and D stands for the diffusion constant. Zero-mean 
thermal noise is characterized by its variance (rj(t,r)r](s,r')) = 2DT5(t — s) S(r — r') where 
T is the bath temperature. The initial conditions are specified in terms of a(r — r 1 ) := 
(0(0, r)0(O, r')) and where we already anticipated spatial translation invariance, hence a(r) = 
a{—r). The potential v = v{t) acts as a Lagrange multiplier. For z = 2 it is easy to see that if 

k(t) :=cxp^-D J^duv(u)^j ~ t F , F = 1 + a - ^ (7) 

and if a = a 1 then eq. (@J is reproduced from Schrodinger- invariance [17]. The Langevin 
equation © may be turned into a field-theory using the Martin-Siggia-Rose formalism. Pro- 
vided that field-theory is Galilei- invariant in the absence of thermal and of initial noise (i.e. 
T = and a(r) = 0) then the two-time autocorrelation function can be expressed in terms 
of noiseless response functions. Precise data on the form of the space-time response func- 
tion in the Glauber-Ising model in 2D and 3D provide strong direct evidence in favour of 
its Galilei- invariance [20] . We concentrate on the case of a fully disordered initial state with 
a(R) = 2ao6(R) where do is a normalization constant. Then it is shown in [17] that 

C(t,s) = a Q JdRR { 3) {t,s,0;R) + DT JdudRR { 3 \t 7 s,u;R) (8) 
4 3) (M, U ;r) := (<f>(t; y)<f>(s; y)4>{u- r + yf)^ - ^^^-U^ (t, s,u;r) (9) 

where the index refers to the noiseless part of the field-theory. Here, the field 4> has the 
scaling dimension x = 1 + a and <f> 2 is a composite field with scaling dimension 2x2 (only for 

frt\ 

free fields X2 = x). The well-known three-point response function 1Z for v(t) = is given 
by Schrodinger-invariance [31] 



7?.g 3 ^(£, s, u; r) = 1Zq (t, s, u) exp 



M t + s - 2u 2 
~2 (s - u)(t - u) V 



q, [ r 2 

2(t - u)(s - u) 



?4 3) {t, s ,u) = e(t - u)Q( s - u) (t - u y* 2 (s - u y* 2 (t - S y x+Z2 (10) 

where ^ = ^(p) is an arbitrary scaling function and M. = \ j(2D) is a non-universal constant. 
Eqs. I|8I9I10|) are the foundation of our analysis of the autocorrelation function. 

Comparing H8I9I10|) with the scaling form we have X2 — x = d/2 — Ac and [17] 



C(t,s)=a y^/ 2 (y-iy Xc <l>(y±±) , dR exp 



Mw r2 



^ (R 2 ) (11) 



where y = t/s. The second term in © merely gives a finite-time correction and may be 
dropped, in agreement with T being irrelevant for T <T C [1]. The form of fc{y) still depends 
on the unknown function <I>(u>). A simple heuristic way to fix its form is to argue that the 
noiseless response function R^p {t, s, 0; r) which describes a response of the autocorrelation 
C(t,s) — (cf)(t)4>(s)) should be non-singular at t = s. This leads to $(w) ~ $o w_Ac as 
w — > oo. If this were valid for all vo, we would obtain the following simple form [17] 

C(t, S )*a <S> ((y+l) 2 / y y Xc/2 (12) 
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Table I - Parameters of the autocorrelation function of the 2D Glauber-Ising model. 



T 


Coo C rc f 


ABE 


0.0 
1.5 


1.65(3) 0.605 3.5 
1.69(3) 0.48 5.5 


-0.601 3.94 0.517 
-5.41 18.4 1.24 



which at least gives the correct asymptotic behaviour as y — ► oo. It has been checked that 
this form is exact for systems described by an underlying free-field theory [17]. 

We now outline how to find the scaling function $(w) in (|ll|l more systematically. To 
achieve this by a dynamical symmetry argument, an extension of the Schrodinger group used 
so far as dynamical group has to be found. Indeed, when considering the dynamical symmetries 
of the free Schrodinger equation (2Mdt — d 2 )4> = 0, it is possible to consider also the 'mass' M 
as a dynamical variable [32]. Then the dynamical symmetry group extends to the conformal 
group in d + 2 dimensions [29,33]. We postulate that, at zero temperature, this conformal 
symmetry is a dynamical symmetry of phase-ordering. Now a quasiprimary field depends on 
three variables <fr = cf>(M,t,r). For conformal invariance, it is sufficient that <fi transforms 
covariantly under the extra generators [29] (for simplicity we also set d = 1) 

d 2 d d d 

V- = -iTTTTTT + ir— , 7V = -t— - 1 - (13) 
dMdr dt dt dM 

We look for the general form of a three-point function — {4> a (f>b<t l c) of a theory without 
noise and which is conformally invariant. Recall that in the MSR formalism, response functions 
are written as correlators of <fi and <j>. The response field (f> conjugate to the field <j> nas the 
'mass' M = -M < [29] and for TZ { Q 3) the 'mass' conservation M a + M b + M c = holds [31]. 
The response function is given, because of Schrodinger-invariance, by eq. I|10[) . Since now the 
'masses' are also considered as variables, we expect $ = ^(p, Ai a , Mb)- Covariance under 
the conformal group means V^TZ^ = N a TZ^ = and leads to 

d d d 2 \ / d Kd d d 2 \ 

A c - a - - + MaTTTT- + E ETl * = 0, A c - a - - + -M hT — — + - — * = 



2 dM a dpdMa) V 2 oM b dpdM 

Xc-2a-^-p—+M a ——+M b —-)^ = (14) 
2 dp dMa dM b J 

hence * = p Xc - 2a - d l 2 K(r)X) with r\ = (M a + M b )p/2 and C = [M a - M b )p/2. We need 
the response of the autocorrelator {<j> a <t>a), thus M a = Mb — M, hence C — 0. Then 

[(2A C - d - 2a) + (A c + 1 - 2a - d/2) 8 V + r]d v + rjd 2 ] K(rj, 0) = (15) 

and we finally obtain the required scaling function (V>o,i are arbitrary constants) 

V(p,M,M) = ipo p Xc '- 2a - d/2 1 F 1 (2A C - d - 2a, X c + 1 - 2a - d/2; -Mp) 

M 2a+d/2 ~ Xc iF 1 (Ac - d/2, 1 + 2a + d/2- \ c ; -Mp) (16) 

Before we can insert this into (|llfl . we should consider the conditions required such that 
the derivation of i|16fl is valid. In particular, it is based on dynamical scaling and we recall the 
condition i-s> £micro f° r its validity (similar difficulties have been encountered before for 
integrated response functions, see [10,11,34]). From 1)100. this means that for small arguments 
p — > the form of the function ^(p) is not given by local scale-invariance. Rather, for p <C 1 
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we expect that the response of the two-time autocorrelation function C{t, s) = C(s, t) should 
be symmetric and especially non-singular in the limit t — s — > [17]. This suggests that 
^>{p) c± ^Qp Xc ~ d / 2 if p < e and ^(p) is given by i|16|) only if p > e where e sets the scale which 
separates the two regimes. The constant Vto is determined from the condition that ^S(p) is 
continuous at p = e. A straightforward but slightly lengthy calculation gives 

$0) = B \{T(d/2) 2 Fi (A c - d/2, d/2; 1 + 2a + d/2 - X c ; -1/w) - j(d/2, Ew)) w~ d/2 
+A [r(A c - 2a) 2 Fi (2A C -d-2a,\ c - 2a; X c + 1 - 2a - d/2; -l/w)w 2a - Xc 

-7(A C - 2a, £w)w 2 "- Ac ] + AE- 2a itr Ac 7 (Ac, Ew) + BE d/2 ~ Xc w~ Xc 1 {\ c , Ew) 

+AE i-2a 2X c -d-2a Xc ^ Ew) 2a-x^ Xc + 1 _ 2 a, Ew) - 7 (Ac, Ew)] (17) 
Ac + 1 — 2a — a/2 

where 7(0, z) is an incomplete gamma- function, E = Me and A, B are constants related to 
V>o,i- Eqs. (|1 111 7|) together give the autocorrelation function C(t, s). This is our main result. 

Some simple consistency checks are easy to perform. First, for free fields Ac = d/2 and 
eq. <|12l) is recovered for A = 0. Second, we find $(iy) ~ w~ Xc for w — > 00 as expected and if 
A ^ 0, we obtain the additional constraint 2a < 1. 

For a non-trivial test, we consider the phase-ordering kinetics of the 2D Glauber-Ising 
model, which we realize through a standard heat-bath rule, and lattices up to 800 2 . We 
consider quenches to temperatures T = and T = 1.5, both in the ordered phase. For 
T = (T = 1.5) we went up to y = t/s = 100 (y = 60), with s = 1600 being the longest 
waiting time studied, and averaged over typically 500 independent runs for the largest lattices. 
While the exponent a = 1/2 is known [3,10], we repeated the determination of Ac and find 
Ac = 1.25(1), in agreement with earlier results [6,26]. Next, we determined the amplitude 
Coo from C(t,s) ~ C 00 (t/s) _Ac '/ 2 as t/s — > 00 which produces a first constraint for the fit 
of the constants A,B,E. A second constraint follows from the observation that at a special 
value y r ef the curves for different values of s cross. We write C(t/ re f) = C ro f. The results are 
listed in table [IJ together with the values of the parameters A, B, E into which we absorbed 
the normalization constant ao- In figure ^ we finally compare our Monte Carlo data obtained 
for T = with several theoretical predictions. Clearly, the waiting times considered are large 
enough to be inside the dynamical scaling-regime. 

Considering first a large range of values of t/s (see figure^,) we observe that although 
the prediction JSJ [22-24] is quite close to the data for t/s small (even so it lies systematically 
above the numerical data), there is a strong deviation for t/s > 3, see [26]. On the other hand, 
the simple approximation (|12f) works well for t/s large but not surprisingly fails for t > s since 
the model at hand is not described by a free field. Finally, local scale-invariance (LSI) as given 
by eqs. (|11I17|) and the parameters of table ^produces a nice overall agreement with the data, 
up to the smallish region t/s < 2. That region is examined closer in figure We remark that 
for t ~ s dynamical scaling no longer holds true [34] (see the inset in figure^)) and we cannot 
hope to be able to find C(t,s) from a dynamical symmetry argument. The approximate 
analytical theories proposed in [13, 35] are only qualitatively correct which indicates that 
OJK-style approximations might not capture fully the quantitative aspects of phase-ordering. 

In figure|3 we present C(t, s) in a more traditional way usually preferred by experimental- 
ists, for both T = and T = 1.5, which makes the simultaneous dependence of C(i, s) on t — s 
and on s explicit. Again, we find a nice agreement between LSI and the numerical data, but 
with larger finite-time corrections to scaling for T = 1.5 than for T = 0. This finding is strong 
evidence that the extension of dynamical scaling to Schrodinger-invariance and further to 
conformal invariance involving also the 'masses' as variables is indeed a true dynamical sym- 
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Fig. 1 - Scaling of the autocorrelation function C(t, s) of the 2D Glauber-Ising model at T = 0. The 
curves are as follows: BPT corresponds to @, app to H'2l and LSI to I|llll7fl . Error bars are much 
smaller than the symbol sizes. The dash-dotted line in a) gives a 2 nd -order perturbative correction 
of @ [13] and the dash-dotted line in b) is the closed-form approximation of [35]. The inset shows 
C(t, s) for s = 200, 400, 800 and 1600 (from bottom to top), as well as the BPT line eq. ©. 



mctry of phase-ordering kinetics for all temperatures T < T c . We recall that the LSI-prediction 
(|11I17|) depends on (f> and <fi 2 being quasiprimary under Schrodingcr/conformal transforma- 
tions [18]. It has turned out that the magnetic order-parameter of the Glauber-Ising model is 
indeed quasiprimary. For the XY-model, however, the spin magnetization S(t, r) cannot be 
identified with a quasiprimary field but the spin-wave approximation suggests that the phase 
variable <f>(t,r) should take that role [17]. 

Summarizing, we have proposed to extend the usual dynamical scaling found in phase- 



1.0 



O 



0.1 



Trm| — i i iiiiiii — i 1 1 mini 

'^':-^ o (a) i 



± « ■ A 



-- BPT 

LSI 
os=1600 
■ s=800 
o s=400 
i s=200 
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10 1 10 2 10 3 10 4 10 5 



0.1 



I I I lllll| 1 I I lllll| 1 I I lllll| 1 I I Mill 

(b) 




10 1 10 2 10 3 10 4 10 5 



t-s 



t-s 



Fig. 2 - Autocorrelator of the 2D Glauber-Ising model at temperatures (a) T — and (b) T = 1.5. 
LSI is the prediction eqs. 1111171 and for a = 1600 the curve BPT eq. is also shown. 
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ordering kinetics to a local scale-invariance by (i) requiring Galilei-invariance at T — and (ii) 
considering the dimensionful 'masses' of the order-parameter and response fields as further 
variables. This has led us to postulate a new kind of time-dependent conformal invariance in 
phase-ordering kinetics. We have derived the explicit prediction eqs. (|1 lllTfl for the two-time 
autocorrelation function. This expression is in agreement with numerical results of the 2D 
Glauber-Ising model and also agrees with several exactly solvable systems [17]. 

* * * 
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